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ABSTRACT 


Phenneacutranstfer-characteristies of a rocket motor 
storage container system have been investigated using 
analytical and experimental techniques. Analytically, both 
closed form and numerical solutions have been developed. 
These solutions may be used to determine maximum tempera- 
tures and temperature gradients within the rocket motor. 
Comparison between theoretical and experimental values of 
temperature are within the estimated experimental uncer- 
Eeamties of GPN3° RF. 9°Lt is proposed ‘that the theoretical 
solutions can be used to thermally optimize container 
design. 

A secondary investigation was carried out to determine 
Biesteaci bility tof wising chelesteric liquid crystals, a 
temperature sensitive material, to thermally map the surface 
oe Che container. |The crystals appear to remain stable 
under desert type conditions and produce brilliantly colored 


displays: of the |itemperature field., 
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IL. INTRODUCTION 


The purpose of this investigation was to develop a heat 
Eransfer model that will allow prediction of the temperature 
fester patton an ‘a c<ontaimer stored rocket* motor placed in a 
hostile thermal environment such as the desert. It is pro- 
posed that such a model would be a useful tool for thermally 
optimizing future container designs. As extreme variations 
in the rocket motor temperature may lead to large thermal 
Stresses in the propellant which could result in fracture, 
or otherwise degrade the performance of the motor, a major 
objective of this study was to design a model that could 
reliably predict the thermal gradient in the motor. The 
predictions would be based on the surface temperature dis- 
tribution, the thermal properties and the geometrical de- 
tails of the system. The model Lay also be used to predict 
a critical temperature range over which the propellant must 
be chemically stable when in a storage situation. The 
upper Limit of this temperature range is referred to as the 
design temperature of the system. As the design temperature 
for most weapon development projects is derived from dump 
storage dowatetens/ a dump storage situation was used to ob- 
GainSthe experimental “data ‘for this*project. 

Several approaches were taken to predict the rocket 
motor temperature distribution from a knowledge of only the 
surface temperature distribution of the storage container 


and the thermal properties and geometrical details of the 
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experimental model. The experimental model used in this 

test was a ian Navy antisubmarine rocket (ASROC) 
motor, filled with dry desert blow sand to simulate the pro- 
peelant, and placed in its. storage container. This container 
system was placed in a dump storage site at the Naval Weapons 
Genter, China Lake, California to simulate a desert environ- 
ment. 

The method of complex temperatures [Ref. 1] was used to 
develop an analytical prediction of the transient tempera- 
Pine eaiteld Ehat exists im a container stored rocket motor. 
The analytical model assumes that heat is transferred only 
tiene radial direction and that the container surface tem- 
DeGrature Variation is sinusoidal with time. Comparison 
between theory and experiment is within experimental uncer- 
tainty when temperature is interpreted as "bulk" temperature. 
ime analytical model is especially useful for studying 
geometrical and thermal physical property effects on rocket 
motor temperature. Such parameter studies have been carried 
Cirpvand the results are presented in a form that will be 
peerul from a container design point of view. 

TRUME WN ReLs..2.and S),.a- computer program for transient 
and steady-state temperature distributions in multidimensional 
Systems, was used to obtain detailed information about the 
Eneamal State of Ehesrocket motor. TRUMP: allows actual con- 
tainer surface temperature distributions to be used as well 
as sinusoidal variations. In addition, both one dimensional 


(radial) and two dimensional (radial and circumferential) 
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heat) transfer were modeled with TRUMP, using both the sinu- 
sondal and actual temperature distributions. The actual 
temperature distributions were obtained from the experimental 
data of the motor container system. 

Comparisons between the experimental values and those 
predicted by the models were in good agreement ,with those 
predicted by TRUMP using the actual temperature distribution 
as the boundary condition being the closest. However, the 
Banwsoidal variations used in.both the analytical moded and 
the TRUMP model are also suitable for design purposes. 

Another aspect of this project was to obtain the storage 
container surface temperature distribution using cholesteric 
liquid crystals, a material that undergoes brilliant changes 
in color over known, well defined temperature ranges. Color 
Slides and movies were taken of the liquid crystals demon- 
Stmhating the feasibidity of using them for on site tempera- 


ture measurements. 
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IL. BACKGROUND 


In 1959 the Naval Weapons Center, China Lake recognized 
Ebeenmeed for a concerted attack on the problem of thermal 
erercerta assignment for new weapon systems. In 1963 a task 
force was established to study the complete environmental 
Cearteria determination problem. The key to this problem 
seemed to be the thermal area in the storage and transporta- 
tion events of any item. It was realized that transportation 
Was a short term situation compared to the storage situation. 
Therefore, the major portion of the life of an item must be 
in storage. There are three types of storage; covered, igloo 
and dump. The dump storage situation leads to the more 
extreme thermal exposure situations which then leads to the 
design temperature. 

As data was not available for the dump storage situation, 
instrumented storage dumps were created at representative 
Pltages won a worbawide basis So that statistical datascould 
Demucelved jon va .variery of ordnance. The first site was at 
China are California, in the middle of the Mojave Desert. 
tieaesete now Was Ene capability to return about 250 channels 
of information on a continuous time-temperature basis (Fig- 
lucent). Obher aretie aud tropical sites.were set up to. study 
extreme conditions. 

The dump storage situation was reproduced to. study the 
extreme Situation. The ordnance was exposed singly, directly 


Sipuaced won the ground, with the long axis aligned in. the 
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north-south direction to allow maximum normal exposure of the 
Pompainer surface to the suns rays. In actual, practice, 
Srduance..is usually stacked and .oriented in other than a 
north-south nee clea eile obey avoiding the extreme situa- 
tion. Ordnance sitting on the ground receives reflected 
madiation from the ground, cannot,.quickly give off heat by 
fConauction, »to the ~ssoil, and is not as apt to be cooled by 

the prevailing breeze; therefore, extreme temperatures 
‘result. 

The most important source of heat to the ordnance is the 
aerect rradiation fromthe «sun; with reflected radiation of 
secondary importance. \For extreme conditions to eccur the 
Windemust be calm (less than 5 knots), the sky clear, and 
the outside air temperature high. After sunrise, the ord- 
Mance skin temperature rises much more rapidly than the 
ambient air temperature; therefore, the surrounding air 
eoolshthe.ordnance;, rather than heats it. 

The rocket motors used for the tests were military sur- 
plus. Even though the material had served its intended in- 
Fleet purpose, it was still representative of new hardware, 
when viewed in a thermodynamic context. When inert rocket 
motors were available, they were used intact; however, in 
most cases, once-fired hardware was used. Thoroughly dried 
desert blown sand, being similar in thermal. properties to 
most propellants, was used to backfill empty rocket motors. 
It was assumed that the thermal response of the sand filled 
motors was essentially the same as actual propellant filled 


motors. 
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LiL. EXPERIMENTAL PROCEDURE 


Although Naval Weapons Center, China Lake had accumulated 
Vast amounts of data in the past, it was decided to instru- 
ment a rocket motor storage container system especially Poe 
tres project. This would allow base data to be taken.exactly 
where it was required. It also allowed variations in the 
System without interfering with one of China Lake's ongoing 
Progects. An ASROC system was chosen for this study. The 
Ourer storage container was 75 inches long with an inner 
diameter of 18 inches and a wall thickness of 1/16 inch. The 
rocket motor was 57 inches long with an outside diameter of 
i_einehes and, a wall thickness of 1/4 inch. Both the con- 
tainer and motor were made of steel, 

The rocket motor storage container system was instrumented 
with 20 gage copper-constantan insulated thermocouple wire 
wWiteh haswan ISA Calibration of + 1-1/2°F over the range -7/5 
to +200°F. Twenty-one thermocouples were originally placed 
on the system with positions indicated im Figures 2 and 3. 
The ambient air temperature was measured with thermocouple 
number 19 which was located in a Stevenson shelter about 60 
feet agi from the system (Figure 4). The thermocouples were 
mounted intrinsically on the moter and.storage, container. 
slwo small holes were drilled.approximately 1/8 inch apart in 
the metal and the individual wires were inserted in the holes. 
The metal was then hammered around the wires until a snug 


fit was obtained. Bead thermocouples were mounted at the 
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Figure 2. Thermocouple Locations on Experimental Systen. 


Five thermocouples were located under the section painted 
Vien Ene liquid exystals., Their flocations corresponding to 
the ones shown above are: #14= #1, #15= #2, #16= #8, #17= #9, 
and #18= #3 (See Figure 3). 
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eenter of the motor and in the air gap. The thermocouples 
located at the center of the motor were supported by small 
pieces of wood several inches from the bead. The use of 
these supports was necessary to keep the thermocouples in 
position when the motor was being filled with sand. After 
all the thermocouples on the rocket motor were in place, 
the rocket motor was filled with dry desert blown sand. 

The wires from the two thermocouples located in the center 
Gf the motor were led cut a’ hole in the end cap. To avoid 
Beetiine of the sand after the motor was in place on the 
site, with a resulting air gap being formed between the 
sand and the motor skin, the sand was compacted by striking 
the sides of the motor with small sledge hammers and then 
Zading “additional sand through the hole in the end cap. 
This was continued until the sand was tightly packed. The 
hole in the end cap was then sealed. The rocket motor was 
Caretally placed in dts storage container (Figure 5) which 
had previcusly been instrumented with thermocouples. The 
thermocouples in the air gap were mounted es at faxine sthe 
Jeszd wire to ther rocket motor. at the desired position and 
then putting a 90 degree bend in. the wire so that it placed 
the bead of the thermocouple approximately 1.5 inches into 
the air gap. Neither the thermocouples in the center of 
the motor nor those in the air gap could be considered accu- 
rately positioned; however, every effort was made to mini- 
mize positioning exrors.. Ali, thermocouple wires were 


located inside the storage container and-were led through a 
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hole in one end. This hole was then sealed. The two halves 
ot the storage container were then bolted shut. 

The outer surface of the rocket motor and the inner and 
outer surfaces of the storage container were all painted 
various shades of haze gray. Weathering had caused the 
ButmEcd surfaces to .uappear fairly rough. This is typical 
Oi the conditions of a-storage dump. From the condition of 
the surfaces, it was estimated that the emissivity was 
approximately 0.9. 

Exnioxr,to Loading the rocket motor into the sterage con- 
tainer, it was decided to apply liquid crystals (See Appen- 
cia A) to part of the storage container surface in order to 
obtain a thermal mapping of the surface temperature at any 
instant of time. Liquid crystals are temperature sensitive 
materials that produce immediate thermal images in a pattern 
of colors which respond rapidly to minute changes in sub- 
strate surface temperatures. A second reason for applying 
the, crystals to the container surface was to determine the 
feasibility of using the crystals under adverse environmen- 
tal conditions (desert atmosphere). Prior to applying the 
erystals, a lo ineh strip of .the storage container, 20 inches 
from one end, was sprayed with two coats of Testors Spray 
Blastnamet No. 12495 Flat Black as,a background, forthe 
eryctalc. sA,one inchs strip of 1.1 different ranges of crys- 
tad, with approximately 1/2 inch. ef black paimt between 
them, was. applied over the black paint. Two coats of each 
crystal were applied, using a small paint brush. The first 


coat was allowed to dry completely before the second coat 
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Maeeapplied, After the erystails were dry, two coats of Rez 
polyurethane (gloss clear plastic coating, interior-exterior 
77-5) coating were applied by brush completely covering the 
Grystals and black painted area. The polyurethane coating 
Was applied to protect the crystals from wind blown sand and 
from the ultraviolet rays of the sun. Ten of the eleven 
crystals had been previously calibrated [Ref. 4]. Using the 
constant temperature bath procedure recommended in Ref. 4, 
R-27 was calibrated and the complete calibration results are 
Shown in Tabiie L. 

The rocket motor storage container system was then moved 
to the China Lake dump storage site. The system was aligned 
in a north-south direction, well away from the influence of 
other ordnance (Figure 6). The thermocouple leads were 
connected through a junction box and underground cable to a 
Honeywell Electronik 25 Recorder which had been calibrated 
to read the thermocouple output directly in degrees Fahren- 
Detiberan accuracy of + 11°F. Thetrecorder was located in 
an air-conditioned building about 60 feet from the system. 

initral data indicated that “the number 7 thermocouple 
was not responding properly and therefore this data was 
Adie ticds Initial color photographs were-taken of the 
id@urd crystals and it was immediately apparent that good 
thermal mappings could be obtained if the crystals were 
Stable under the adverse desert environment. The brilliance 
Gf the colors ‘exhibited by the crystals under the bright 


desert sun was much better than had been expected- The 
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PABGE 1 
Calibration of Liquid Crystals 
Calibration: Bath 


NCR Manufacturer's 2 Goats 
Desig. Color Change Responses Liquid Crystals 
She 2G 

R-27 

Red IF eg) 2 O-tze 

Green ZE00 28.04.55 

Blue 30.0 28. 1a 
R-33 

Red Br). 0 32 FED 

Green 3 ty 16 33 ico akiD 

Blue 361.0 34.2.5 
R-37 

Red ey sO) 36.24. 5 

Green 38.6 37.1+.5 

Blue 40.0 38.0+.5 
R-41 

Red SG) 40.3%.5 

Green 42),0 41.0% 55 

Blue 44.0 42 04.5 
R-45 

Red 45,0 42.8+.5 

Green 46.6 43.6+.5 

Blue 48.0 44.3+.5 
R-49 

Red 499°.10 AO Ito 

Green 50.6 LeTiw dite 5 

Blue 5272.0 48.44+.5 
R-53 

Red 53.0 50). 5335 

Green 510 D2) ste 

Blue ; S670 33.3.5 
R-56 

Red 316.50 I398 45 

Green DG DIO" Os <5 

Blue 59.0 DOg oa sD 
R-59 . 

Red 59" .0 BIG Stee 

Green 60.6 SWAPS aera) 

Blue 62.0 D869. 5 
S-62 

Red 62.0 607. Ia. 5 

Creen, 62.6 60.44.5 

Blue 63°10 60.94.5 
S-64 

Red 64.0 60.9.5 

Green 64.6 61.44+.5 

Blue 65.0 627755 
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system was allowed two weeks to reach a periodic steady 
meaee before additional photographic data was obtained. 

[=eeusave photographic data was collected on 27 and 28 
July 1972 after two weeks of exposure to the desert environ- 
ant. Both super 8 mm and 16 mm color movies and 35 mm color 
Berdges were taken of the liquid crystals. -No colored filters 
were used on any of the cameras, although standard haze fil- 
ters were used to take the super 8 mm movies and most of 
the 35 mm slides. 

Atethisntimes,-a»second storage, container, piewene with- 
out a rocket motor inside, was instrumented with intrinsic 
thermocouples in the same manner as the previous container. 
As only three data channels remained open on the recorder, 
only three thermocouples were applied to this new container. 
The three thermocouples were applied at the 0300, 0900, and 
1200 positions at the midpoint of the container. Ehicmcon= 
tainer was set end to end with the system that was already 
anepilace at the site. The purpose of this study. was to 
Soeermaine 1f the inelusion of:the rocket, motor in, the con= 
Painer had a significant effect, on, the. surface, temperature 
of the container. Thermocouple #7 was connected at the 0900 
Pesttton, £23,.at the 1200 position, and #24, at the 0300 posi- 
tion. It was immediately apparent that thermocouple #7 was 
Continuing to give unreliable readitigs and therefore the 


data taken on channel #7 was neglected. 
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Vicor THEORETICAL ANALY S LS 


As ONE-DIMENSIONAL ANALYTICAL MODEL 
the first step was to try to devise an analytical model 
that Peaia simulate the actual rocket motor storage contain- 
ereexperimental system. The first simplifying assumption 
was that the storage container temperature could be modeled 
by a sine wave which had a period of 24 hours. A compari- 
son of the sinusoidal variation to the average (bulk) stor- 
age 2onte ane temperature [obtained by averaging the four 
thermocouple readings on the surface of the container (l, 
Swe. and 12) as shown in Appendix.D] is given in Figure 7. 
The method of complex temperature as presented by Arpaci 
[Ref. 1] was used to find the steady periodic solution of a 
body experiencing a periodic sinusoidal disturbance. A com- 
plete analytical derivation is given in Appendix B. The 
general heat conduction equation in cylindrical coordinates 
was sthe basis for this derdvation. It was assumed that 
there was one dimensional radial heat flow with no conduction 
imeene axial or circumferential directions, that no sheat 
sourees existed in the model, that the rocket motor storage 
container system was infinitely long, and that the sinusoid- 
ally varying surface temperature was spacially uniform over 
the entire container surface. The storage container tempera- 
ture is assumed to vary as 
a (Ty - T,) Salint “iyi en a) gt 
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remhergeea@ (8) 


where Ty = maximum bulk temperature of the storage 
container 


Ty = average bulk temperature of the storage 
container 

W = frequency of the sinusoidal variation 
(27/24 hours) 

t = time (hour's) 


It was assumed that all the thermal properties remained 
constant over the temperature range of the problem. The 
effective heat transfer coefficient, h, across the air gap 
between the storage container and the rocket motor combines 
the heat transfer effects of radiation, convection, and con- 
duction into one coefficient. The radiation coefficient was 
linearized by assuming constant representative temperatures 
in the equation 
te aye te) 


a 2 


where 0 is the Stefan-Boltzmann constant and es is the 
Fadiation exchange factor. The .convection coefficient is 
the effective conductivity aes obtained from the peeke 
Hann correlations [Ref. 5],.divided by the width.of, the air 
rape cil the analytical model, the, effective conductivity 
was assumed to equal the conductivity, thereby treating it 


as pure conduction and giving the equation 


Baer ne Heo 


An Anitial eendition was not specified in this deriva-= 


tion as the only concern was the steady-state, periodic be- 


havior. The steady-state solution is (Appendix B) 


UGE) plates Es (a EBEa 7(ae) 
Oe, &) = 5 a = ae sin(Wt+6*) 
= SoG ee 
M A Soe bake 
MiGeety = 6 Siu (or + oF) 
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where T(r,t) Ehe temperature of a point xr in the rocket 


Moto ae Eine 


Wr 
a= i = conduction parameter 
E Bon aU 
E = ae dimensionless distance from the center 
© of the rocket motor 
rhe inner radius of the rocket motor 


r= distance from the center of the rocket motor 
k : as 

a me = thermal diffusivity 
Pp = density 

k = thermal conductivity 

c 


= specific heat 


BER = real Bessel Function 
BEi = imaginary Bessel Function 
ee ore te ER ape —— BT Ca) 
R fe) 1 1 
28 
YenHibes Ga) esr BEA ie) = —-_ BERe (a) 
a ; fe) 28 af: /2p IL 
hry 
B = gy Biot modulus 


Bid (ab)x, = BER (ae)x. 
BER | (a&)X, + BEi | (a&)X, 


6* = Pane 
Two computer studies were done based on the steady state 
solution. The first study was a completely dimensionless 
Situation which served as a parameter study of the effects 
©f varying a and 6 on the temperature and the time lag of 


the temperature at various positions in the model. 


wr 
a = A = conduction parameter 
and a 
hr) 
B = x = Biot modulus 


Paraneter! a was’ varied from 41.0 to 5.0) and — was’ varied 


from 0.1 to 100. These were the only values studied, as 
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only values within this range are of interest in this type 
problem. The computer program and its output are given at 
the end of Appendix B. The output lists the following 
values: 

1) a, the conduction parameter 

eG, thes Brot. modulus 


3) —&, the non-dimensional distance from the center of 
the motor 


4) 6*,the time delay between the maximum storage con- 
tainer temperature and the maximum temperature 
reached at..the point) of (interest in the motor 


5) 0 .,the relative amplitude of the maximum tempera- 
ture at the point of interest compared to the 
maximum temperature of the storage container 


The time delay is given in radians, where 27 radians equals 


one complete cycle. A graph of the time delay versus § for 


a@ecomstant value of "a" is given in Figure 8 at three dif- 


ferent positions within the motor. A graph of the relative 


amplitudes of the temperatures versus §$ for a constant value 


Sea 15 givyenin Figure 9. Lt was meted that -for a con= 


stant value of "a", the time delay decreased as £8 became 


tareer. As the point of interest approaches the center of 
the rocket motor, the time delay increases. The relative 
amplitude of the temperatures also becomes larger as 8 is 


Mmereascdawhen the value of Va" is held eonstant. Pf 6 is 


held constant and "a" is varied, the time delay increases 


and the relative amplitude decreases as "a" increases. 


The second study was obtaining the analytical solution 


to the particular rocket motor storage container system 
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studied at China Lake. The thermodynamic properties of dry 


sand were obtained from Rew. .Oaas 


p= 194 aGwileuseee 
f= 0.188 BIU/bx ft °F 
ees ORT) tbr °F 


Substituting these values and using 1440 minutes (24 hours) 
@s aycomplete cycle, the parameters a and §$ were calcu- 


lated For this model as 


ie, 
a= =2543 
a 


where rea 5S./5 anehes, the inner radius of the rocket motor. 
hr | 
Bp = aie 2090 
where h = hoon + Apap 
k 
and hooy = Got = 6-48 x 197* TU 
': pre et i 
where) -Ar.='2.94-inehes,. the distanee across: the air gap 
-2 BTU 
2 2 BTU 
ee sooth er an) 1.09 4520 
RAD 1-2 i. 2 il 2 pyspedosp 


where oO is the Stefan-Boltzmann constant, uae, is the 


Fadiation exchange factor which for this geometry is 


= 0684 
set) 


when Eye eg = aos eo = 165.0; La ke 8.94 
therefore h = 1.15 BTU/hr f£t2°F 
The average surface temperature of the storage container 


was tound to be PO4°F for a particular day at China Lake, 
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with a maximum temperature of 138°F. These values were ob- 
tained by averaging the readings of thermocouples 1, 8, 10, 
and 12 as shown in Appendix D which give the bulk tempera~ 
EW ee), | 

The temperatures of seven positions within the rocket 
motor were calculated and the results are printed at 30 
Minute intervals for one complete cycle in Appendix B. A 
graph of temperature versus time was plotted by the computer 
Showing the relationship between the surface temperature of 
the storage container (TINF), the temperature on the outer 
Skin of “the rocket motor (TEDG), and the temperature at the 


eemeer of the motor (TCEN).»§ This graph is Figure 10. 


B. TRUMP MODEL 

ihe recket moter storage container system at China Lake 
Was modeled on TRUMP, a numerical conduction code, (See 
Appendix C for a description of the TRUMP program) to pre- 
dict the temperature at any point in the system from a know- 
ledge of the storage container surface temperature variation, 
the thermal properties and the geometrical details of the 
System. Two models were used to simulate the rocket motor 
storage container system and several variations of each 
model were investigated. 

The first model assumed one dimensional heat transfer 
(radial). The system was modeled as two infinitely long 
Concentric steel cylinders, the inner of which was filled 
with dey sand. A’ 2.94 inch air gap separated the cylinders. 


The model was subdivided into concentric volumetric elements 
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with representative nodal points as given in Figure 27, 
appendix C. It was assumed that the storage container sur- 
face temperature was spacially uniform. From the data given 
mrmeAppendix D and the observation of the liquid crystals' 
thermal mapping, it was obvious that the temperature distri- 
bution on the storage container was not spacially uniform. 
In order to simulate a spacially uniform condition, the 
readings of the thermocouples located at the 1200, 0300, 
mao, and 0900 positions (#1, 8, 10, and 12) were averaged 
and this average value of the surface temperature (referred 
to as the bulk temperature) was used as the spacially uniform 
temperature distribution. Two methods were used to describe 
the container temperature. The first method used the maxi- 
mum bulk temperature (138°F) and the average bulk tempera- 
ture (104°F) of the storage container to generate a sine 
wave with a period of 24 hours (1440 minutes). The second 
method took the bulk temperature readings at two hour inter- 
Vals and fed this data into the TRUMP program in a tabular 
(temperature versus time) form. The version of TRUMP used 
in this problem was Limited to a table length of 12¢tabular 
values. TRUMP interpolated between the tabular points. 
Figure 11 compares the actual bulk data with the sinusoidal 
approximation and the interpolated tabular values. 

Several assumptions were made to simplify the solution 
of this problem. As the thermocouple data from the storage 
container gave an average value of the temperature across 


tne 1/16 inch steel’ wall, node’ 12 was modeled as a zero 
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volume boundary node with a known temperature eT ee upon 
Meee Lt was also assumed that heat transfer across the air 
gap occurred by radiation and conduction alone. Free convec- 
tion effects were initially neglected. This assumption was 
later modified to investigate the free convection effects. 
All surfaces of the storage container and the outside surface 
of the rocket motor were painted various shades of haze gray 
and it was estimated that the emissivity of these surfaces 


was 0.9. The radiation exchange factor, J _2for this 


1 
model was the same as that for the analytical solution Cae 
= 0.84). It was also assumed that there was perfect thermal 
contact between the rocket motor and the sand that filled 
Sit, This meglects the possibility that the sand might 
comehtiy settle after being on the site for a long period of 
time. 

The second model assumed two dimensional heat transfer 
(radial and circumferential). The same physical model was 
used as in the one-dimensional case with the sole exception 
that 48 nodes were used instead of 12. The representative 
nodal points and an example of the thermal connections from 
one of-the nodal points are shown in Figure 28 in Appendix C. 
The four nodes on the surface of the storage eae eke 
modeled as zero volume boundary nodes. The sinusoidal and 
tabular representations were used to describe the surface 
temperature of the storage container at each boundary node. 
metual data taken at each position, rather than bulk data, 


were used as the input data for these representations. 
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The same assumptions made in the one dimensional case 
were also applicable to the two dimensional model. A com- 
puete discussion,of, the calculation of the radiation ex- 
change factors in the ve dimensional case is given in 
Appendix C. 

mie effect of natural convection was studied in both 
the one and two dimensional models. References 5 and 7 
Bive correlations between the Grashof number based on the 
gap width and the effective thermal conductivity. The 
Grashof number was calculated from the equation 


6 eB(AT) 6° 


Gr = uw? 


where 6 Width of the air gap 


> 
=) 
" 


maximum temperature difference at any instant 
of time in the air gap 


B = + where T = 565°R 


Pie = 565°R, air has the following properties 
a 


i} 


p 0.07 1bm/ft 


U 0.046 lbm/hr-ft 


The maximum Grashof number for this experiment was calcula- 
fea to be 1.25 x Too. The diameter ratio was approximately 
i> and the log Gr = 6.1. From the Beckmann correlation 

k 


[Ref. 5], the effective thermal conductivity ratio Ga was 


approximately 3.2. WUWsing the Liu correlation [Ref. 7] 


ko ages 0.278 
eet? ei serPe is 


where the Prandtl Number = 0.707. An effective thermal con- 


ductivity of 4.0 was assumed as the average value of these 
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Emo correlations and it was used to study the effects of 
free convection. This change was placed into the TRUMP pro- 
gram by increasing the value of the thermal conductivity of 


PeEeby a factor of 4 in each of the TRUMP runs. 
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Ve RESULTS 


aa) ANALYTICAL MODEL 

Using the sinusoidal temperature distribution as an ap- 
proximation to the actual average experimental data as 
shown in Figure 7, comparisons were made between predicted 
temperatures and actual temperatures for two radial loca- 
Eions in the rocket motor. Figure 12 compares the results 
Ba the surface of the rocket motor and Figure 13 does the 
comparison at the center of the rocket motor. An uncertain- 
ty analysis is given in Appendix E which establishes the 
uncertainty bounds for both the predicted and the actual 
temperatures. These uncertainty bounds are included in 
Pioures 12 and 13. 

It is readily seen from Figure 7 that a sine wave was 
Bot an ideal fit as an approximation to the experimental 
taba, as lt varies as much as 20°F during part of the cycle. 
However, it was also noted that the sine wave closely oer 
mated the experimental data during the heating phase of the 
cycle and only during the cooling phase were there large 
variations. As the main purpose of this study was to design 
a model that would be useful in optimizing storage container 
Besaten, the errors “dn. the cooling phase are not critical as 
long as the temperatures reach the same minimum point before 
beginning another cycle. Figure 12 shows that the maximum 
surface temperature of the rocket motor predicted by the 


analytical model is a good approximation to the actual 
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Papertimental data. Again it is noted that, in actuality, 
the motor cools faster than the predicted value. The maxi- 
mum difference in temperatures on the surface of the motor 
moval 5° F.. Figure 13 shows that the predicted value and the 
Saperimental value of the temperature at the center of the 
Bocket Motor were in close approximation except during the 
early stages of the cooling phase where a maximum tempera- 
mine variation of about 5°F occurred. 

One of the reasons the system cools faster than predicted 
Fould be the light breeze that. is usually evident in the 
Carly afternoon hours at China Lake that is not present dur- 
ing the morning. No attempt was made to shield the system 
moom the wind to study the effects of a light breeze on the 
surface. temperature of the storage container. 

Mnother point not taken into account) by the analytical 
modcin bs the fact that, the time delay at any point in the 
System is not constant throughout the day as predicted in 
Haare! LO, but varies as given by the data in Appendix Dis 
Time delays between the peak temperature on the container 
Piigace and the peak temperature at the center of the rocket 
motor vary from about 250 to 400 minutes, whereas the low 
temperature on the surface of the container and the low ten- 
perature at the center of the rocket motor vary from about 
#50 to 250 minutes. . The analytical model predicts a constant 
Variation of 388 minutes at the center of the rocket motor 


and i259 minutes at the surface’. 
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B. TRUMP MODEL 
1. One Dimensional 

Four variations of the one dimensional TRUMP model 
were investigated and compared to the experimental data. 
Peeunes 14 and 15 compare the TRUMP predictions to the act- 
Mal experimental data at the surface and the center of the 
rocket motor,respectively. The TRUMP variation used for 
this comparison modeled the storage container temperature 
W2th tabular data (See Figure 11) and assumed convection 
was present a = 4.0). The uncertainty analysis (Appendix 
E) established the uncertainty bounds for both the experi- 
mental and the analytical data in these Figures. The varia- 
tion between the bulk temperature predicted by TRUMP and 
the experimental data closely matches with only two experi- 
Hemeal points in Figure 14 falling outside the uncertainty 
bounds for this one dimensional model. Figure 11 shows 
that the tabular data that TRUMP interpolates is a good 
approximation to the averaged experimental data. At the 
center of the motor, as shown in Figure 15, all experimental 
points fall within the predicted error bounds. A comparison 
cA,chne Lour one-dimensional IRUMP variations are given in 
Hieures 16 and 17 at the surface and the center of the rocket 
motor respectively. It is clearly seen from these Figures 
that the convection assumption results in an increase of 2°F 
in the maximum temperature and a decrease of 2°F in the 
minimum temperature on the surface of the rocket motor. This 


Eeuperature ehange drops, to,+,.d452F at the center of the 
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motor aS Shown in Figure 17. The differences between the 
sinusoidal approximation and the tabular approximation of 
the experimental data was clearly shown in Figure ll. The 
aaea in Figures 16 and 17 can be easily correlated to that 
in Figure 11, thereby explaining the differences in the pre- 
dicted values. 
2. Two Dimensional 

Four variations of the two dimensional TRUMP model 
were investigated and compared to the experimental data. 
Comparisons of each TRUMP variation to the experimental data 
are given in Figures 18 and 19 for node 8 (located on the 
skin of the rocket motor at the 1200 position) and node 1 
(at the center of the rocket motor) respectively. These 
Figures show that the TRUMP variations that used tabular 
data to model the surface temperature of the storage con- 
tainer predicted temperatures that more closely approximated 
the experimental values than were those predicted by TRUMP 
variations using sinusoidal data to model the surface tem- 
perature. Appendix D shows that all points on the surface 
of the storage container reach their minimum temperature at 
the same time; however, these points reach their maximum 
temperature as much as 200 minutes apart. Whereas, all the 
points on the surface of the storage container are in phase 
at the minimum temperature, they rapidly become out of phase 
as the container temperature rises. This varying phase 
SHitth maces 2: difficult to model the four boundary nodes 


with sinusoidal approximations which must have constant 
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Phase shifts. Sizable errors in the input data during some 

Parts of the cycle were caused by these varying phase shifts. 

miese exrors in the input data led to the variations in the 

predicted temperature values. As noted in the one dimen- 

eronal section, the inclusion of convection effects does 

mot produce large variations in the predicted temperatures. 
Figures 20 and 21 show the actual temperature dis- 

Preabutions onthe surface. of the—storage container and on 

the surface of the rocket motor respectively at maximum 

bulk temperature compared to a two dimensional TRUMP program. 

The TRUMP variation used for this comparison assumed no free 

Comvection in the air gap and used tabular data to approxi- 


Mate the surface temperature of the storage container. 


C. GENERAL 

A comparison was made between surface temperatures on 
the storage container that contained the rocket motor and 
the storage container that was empty. The low temperature 
was about 4°F colder in the empty container, whereas the 
high temperature was about 4°F higher.on the container that 
contained the rocket motor. The empty container had a 
faster’ response time than the one containing the motor. 
The differences in heat capacities, eee eifiects, and 
natural convection all contribute to the changes in tempera- 


ture noted. 
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Pogaure 20% temperature Distribution at Surface of Storage 
Container at Maximum Bulk Temperature at approximately 
"D000 on’ 2-August’ 1972. 
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Prgure 20: Temperature Distribution at Surface of- the 
Rocket Motor at Maximum Bulk Temperature at approxi- 
mately 1500 on 2 August 1972. 
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Pe LEQULD CRYSTALS 

The encapsulated cholesteric liquid crystals applied to 
Ene Surface of the storage container gave brilliant colors 
under the intense desert sun. These colors were much 
clearer and brighter than the same crystals viewed under 
Piperatory Lighting conditions. The Jiquid crystals photo- 
graphed well in both the color movies and the color slides. 
The movies showed by time lapse photography the rapidly 
changing surface temperature of the storage container. Two 
Sample color prints made from the color slides are enclosed 
asenigpures 22 and 23 to show the brilliance of the colors 
and the feasibility of obtaining data from color photos. 
The only photographic problem encountered was the intense 
Beptection of the sunlight from the polyurethane film. This 
problem was partially overcome by taking the photographs 
from angles where the reflection was less intense. Quali — 
tatively the liquid crystals were not adversely affected by 
the sun's rays after two weeks of desert exposure. No 
accurate quantitative determination was attempted; however, 
rough approximations were made at the site. These approxi- 
mations were made by. noting the color exhibited by a crystal 
ata Certain time and then comparing the calibration of the 
crystal (Table 1) to the temperature recorded by the thermo- 
couple Jocated directly beneath the region of color change. 
The readings were within + 2°F, which was very encouraging, 
especially considering the approximations made while taking 


these measurements. Although photos were taken only during 
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the initial two weeks of the study, on site observations in- 
meeate that the crystals are still showing brilliant colors 
meeer 3-1/2 months. Preliminary ere indicates that the 
polyurethane film did protect the crystals from decomposing 
from the sun's rays and from being’ worn away by the wind 
blown sand. 

ft was noted that the surface temperature of the storage 
container under the liquid cyrstals reached temperatures up 
Remo § higher than a similar point not under the crystals. 
This 15°F difference was only evident when maximum tempera- 
tures were obtained. During sunlight hours the temperature 
under the liquid crystals was always somewhat higher; how- 
ewer, at night both temperatures were equal. The difference 
in the container surface temperatures led to a difference of 
ma fvon the surface of the rocket motor and 1°F at the center 
Or the motor. Lt is believed that the difference in emissi- 
Wttes of the gray and black surfaces resulted in the differ- 


€nce in container surface temperature. 
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Vie CONCLUSIONS 


Prom the LestLtes Of Ehis investigation, the following 
conclusions were drawn: 

fm xtehouch a sine wave is"not a perfect fit for the 
eeperimental data at all points, it is useful in predicting 
bulk temperatures in the rocket motor, especially if only 
Pwe high and low bulk temperatures are of concern. This is 
especially true in the one dimensional case. In the two 
dimensional case, the problem of phase shift variations 
make the method of sinusoidal variation less desirable, 
atenough still useful. 

Zee he cimgiacton Of the actual data by a ‘table’ of 
temperatures gave the most accurate predictions of the 
experimental data. This method should be used whenever tabu- 
lated data are available; however, this will generally not 
be the case for design work, in which case the sinusoidal 
approximations must be used. 

a ane ti lexibilacy of both the analytical and computer 
simulations allow the changing of many parameters. The re- 
sulting effects of these changes on rocket motor temperature 
can be studied with the models, i 

4. The convection assumption for this system resulted 
in only small changes in temperature and can be neglected 
when predicting design temperatures. Either the Liu or 
Beckmann correlation should be used to determine if convec- 


tion can be neglected in a particular system. 
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>». the use of an empty storage container to obtain sur- 
face temperature data is a good approximation to using one 
Mierie a rocket. motor inside. 

Sea suet Coe Dien tOn use, igquid,erystals for thermal 
mapping under desert conditions. Color photography with 
Standard equipment gives excellent results since brilliant 
colors were observed. 

Pe ne weidudd .coyctads appear to, be stable for at Jdeast 
two weeks under the desert conditions when protected with a 
polyurethane coating. 

bow oe. application. of, thee liquid erystal, system, tio, the 
Surface of the storage container resulted in large increases 
in the surface temperature of the container throughout the 
hottest part of the day. Care must be taken in applying 
and interpreting thermal readings from liquid crystals when 


exposed to radiant heating. 
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“VII. RECOMMENDATIONS 


From the results of this basic study, the following 
Pe Pniendattows for future work -are offered: 

ime  bowrerine tie sesultia.of this project, a second 
rocket motor storage container system should be instrumented 
with the following changes: 

a. Liquid "crystals ‘should not ‘be applied to the 
System used as the experimental model. As steel is a good 
Peermal conductor, axial conduction on the surface of the 
Storage container may be significant. Heat flow from the 
area where the crystals are applied may lead to higher than 
normal temperatures at other points on the surface of the 
container. 

b. The rocket motor should be weighed before and 
atter the loading of ‘the dry sand so that'an accurate deter- 
mination of the density of the propellant simulant can be 
determined. 

c. Four additional -thermqcouples should be located 
on the surface of the storage container to better enable 
the averaging of data. At present, the #1 thermocouple 
which was used as the average temperature reading of the 
Eap quarter ofethe (surface ’of theycontainer,:ini actuality is 
its hottest point; likewise the #10 thermocouple was used 
as the average temperature of the bottom quarter of the 
Mevacew tmvactuatity at's the coldest ‘point. For averaging 


data, it is recommended that thermocouples be placed at 0130, 
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0430, 0730, and 1030 and the quarters of the system be di- 
meacd at O300, 0600, 0900, and 1200 to give a more realistic 
bulk temperature. Thermocouples at 1200 and 0600 will pro- 
vide the maximum and minimum temperature of the system. 

2. The TRUMP program should be rerun in both the one 
and two dimensional form, varying the mesh sizes to deter- 
mine the optimum number of nodes. 

3... A long’ term study of the effects of the: desert en- 
Vironment on liquid crystals should be done. fThe crystals 
should be calibrated before being placed in the desert and 
Peent brought,sto,aclaboratory for recalibration at specific 
intervals. 

4. Several modifications should be made to the TRUMP 
program to make it comparable to the version used at Law- 
rence Radiation Laboratory. The variable conductivity sec- 
Elon ~(BLOCK: 2) .andédthe-PLOT subroutine (BLOCK 11) need to 
be corrected. The TIMEP subroutine which allows the setting 
of the problem time interval between data output should be 
added to this version of TRUMP. It would also be advanta- 
eeous to increase the amount of tabular data ithat could be 
read in as boundary temperatures. 

5. From an academic standpoint, the effects of free 
convection in an air gap with varying boundary temperatures 


should be investigated. 
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APPENDIX A 


im-prodwetion to Liquid Crystals 


Liquid crystals were first discovered in 1889 by Reinit- 
zer [Ref. 8] and the investigations of Lehmann which con- 
Ernded to 1915. Liquid crystals were considered to be 
mimoeratory curiosities with no scientific or practical merit 
until the 1950's. They share some of the properties of both 
Prauids and crystals; for example, a typical liquid crystal 
Substance scatters light in symmetrical patterns and re- 
flects different colors depending on the angle from which 
it is viewed. Studies in the last few years have helped to 
clarify the unusual molecular structure of liquid crystals. 
Many applications arise from their ability to detect minute 
fluctuations in temperature, mechanical stress, electromag- 
metic radiation and chemical environment by changes in their 
Color. 

Liquid crystals are divided into three classes; smectic, 
nematic, and cholesteric, depending on the degree of spatial 
arrangement of the molecules in the mass of the material and 
ewe type of the material [Ref. 9]. In this project only 
cholesteric liquid crystals were used and therefore only 
their properties will be memtioned. The molecular structure 
of cchotesterie Iigqguid cexrystals is, characteristic of the 
esters of cholesterol (Figure 24). The molecular layers 
are very thin with the long axis of the molecules parallel 


te the plane of the layers. The individual molecules are 
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Peeure 24: Molecular Structure of Cholesteric Ester 
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Hagure, 25.6 Light Reflection from Liquid Crystals. 
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masteally £lat, with a side chain of methyl groups (-CH,) 
pen,ecting upward from the plane of each molecule. This 
comtiguration causes the direction of the long axis of the 
Merecules in each layer to be displaced slightly from the 
Gerresponding direction in adjacent layers. This displace- 
ment, which averages about fifteen minutes of arc per layer, 
is cumulative through successive layers, and the overall 
displacement traces out a helical path. 

ite molecular structure of cholesteric, liquid crystals 
Baves rise to many peculiar optical properties. If linearly 
peearazed Light is transmitted perpendicularly to the molec- 
wear Jayers, the direction of the electric, vector of the 
Dient will be rotated to the left in a helical path. There- 
Bore, the plane of polarization will also be rotated to the 
Pett, through an angle proportional to the thickness of the 
transmitting material. Liquid crystals are the most opti- 
Cally active substances known. . Another strictly crystalline 
Sptrteal property exhibited by cholesteric liquid, crystals 
fewer renular diehroism. When ordinary white light is ineci- 
@ent to a cholesteric material, the light is separated into 
two components, one with the electric vector rotating clock- 
wise, the other rotating counterclockwise. Depending on 
the material, one of these components is transmitted, and 
Ene other as, reflected. It is this property that gives the 
enolesteric phase ats iridescent color, when it is .illumina- 
bed by white ticht. The particular combination of colors 
depends on the material, the temperature, and the angle of 


Bue incident Light. 
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ihe moLecularcstructure of a cholesteric substance is 
very delicately balanced and is easily upset. Any small 
disturbance that interferes with the weak forces between the 
molecules can produce marked changes in optical properties 
Buen as refbéction, transmission, birefringence, circular 
Seehnroism, optical activity and color. The most striking 
@ptical transformation that occurs in a cholesteric sub- 
Stance, in response to small changes in its environment, is 
ExeeweVariation of colorewith temperature. The crystal lat- 
tice is disrupted by the thermal vibrations giving succes- 
Sive transitions between the solid, the mesophase, and the 
isotropic liquid with rising temperature. The change from 
the three dimensional order of the crystal lattice to the 
disorder of the isotropic liquid occurs via one or more in- 
termediate states, each of which has a particular tempera- 
lmeeerance até whitch iti ds! stable [Ref. 10]. 

A cholesteric liquid crystal system responds to changes 
in temperature by sequentially passing through the complete 
Visual’ Spectrum (red through violet) in. fractions or, multi- 
degrees, depending on which cholesterol esters comprise the 
formulation. This color phenomenon is reversible and has 
been reported to function over a temperature range of -20°C 
EO 250° CG, PAS weay® importants ipoint tormotesds ‘that, at a cer- 
tain temperature a given material or combination of materials 
Wid always Gxhibit’ the same color.: Also, the rate of change 
from color to color as well as the exact temperature at 


which the specific color changes occur are invariable. By 
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mex icholesteric substances. in variows’ proportions, any 
desired temperature combination can be obtained. The thick- 
ness of the cholesteric film does not affect the predominant 
men iength of the reflected light; the light becomes cir- 
Sularly polarized’ [Ref. 11]. 

The colors scattered by ee liquid crystals represent 
Smly a fraction of the incident light (Figure 25). The re- 
Maining portion of the incident light is transmitted by the 
miquid crystals. Therefore, an absorptive black background 
Must be used to prevent reflection of the transmitted light, 
thereby enhancing the resolution of the scattered colors or 
wavelengths reflected by the liquid crystal systen. 

The cholesteric liquid crystal systems often present a 
number of problems due to the fact that they are viscous 
liquids. Some problems associated with the handling and 
the use of these materials are: 

me he tendency of the liquid crystal system to flow 
euring application tan cause variations in applied film 
thickness. This may result in non-uniform thermal patterns. 

2. -Direct exposure of liquid crystals to adverse en- 
Vironmental effects can cause variations in their sensiti- 
Wity and deteriorate their color response in a few days. 

These problems can be partially overcome by using an 
encapsulated liquid crystal material system. The capsules 
are 20-30 microns in diameter and are a water-based slurry 
Suitable for application by conventional coating techniques 


such as brushing or spraying. 
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Encapsulated liquid crystals offer several advantages: 

me they convert the liquid crystal system to eae pseudo- 
Sotid, which provides for easier handling, application, and 
use. 

2. | ithey provide Fonger)shelf life by minimizing surface 
contamination and giving protection from ultraviolet light 
PWetyol2). 

one they sexhibit ‘relatively ounlimited fatigue life. 

4. They reduce the eee dependence of the color 
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APPENDIX B 


Analytical Solution 


The method of complex temperature as presented by Arpaci 
tnen. 1) was used to find the steady periodic solution of a 
body experiencing a periodic sinusoidal disturbance. The 
general heat conduction equation in cylindrical coordinates 
was the basis for this derivation. It was assumed that no 
neat Sources existed in this problem, that the rocket motor 
storage container system was infinitely long, that there was 
mo heat conduction in the axial or circumferential direc- 
tions, and that the container surface temperature was spa- 
etally uniform. Figure 26 gives a basic sketch of the systen. 
The assumptions reduced the heat conduction equation to 


a o(r =) 


r or 


a aw 
Gent (1) 


where r is the radial distance from the center of the rocket 
motor, TIT is thie temperature of the rocket motor at time t 
ama position r, and ad-is the thermalediffusivity, a property 
Soe the conductines material. 


k 
ee rae (2) 


where k is the thermal conductivity of the conducting 
Material, p is the density of the material, and c is the 
Specific heat. All thermal properties were assumed to be 


constant over the temperature range of this problen. 
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Figure 26. Analytical Model of Experimental System 
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Tr, 


The boundary conditions used in this derivation were 


er = 0 at re == 0 
and i = - a (T - LD at r = To 


where ar isthe tuner radius of the xocket motor. 


T,, is the known storage container temperature which is 


assumed to vary as 


t= (T7 T,) sin wt + T, 
where 
Tu = maximum bulk temperature of the storage container 
T, = average bulk temperature of the storage container 
; ; ig A 27 
W = frequency of the sinusoidal variation (sq PERE, 
= = time 


bh as the effective heat transfer coefficient across the 
air gap between the storage container and the rocket motor. 
It combines the heat transfer effects of radiation, convec- 
Evon, and conduction into one coefficient. The radiation 
coefficient was linearized by assuming constant temperatures 
(T, > Ty); representative of the average temperatures expec- 


ted in the problem, in the equation 


2 78 
ane Meet tgh et To) CT, + Ty) 
where o is the Stefan-Boltzmann constant and wae is? the 
Fadiation exchange factor between surfaces 1 and 2. The 


convection coefficient is 


ee SS 
Hoon mS 


where ko ie the effective), conductivity of air as obtained 
from the Beckmann and Liu correlations [Ref. 5 and 7] and 


Deteoutpe width of the airgap. In the analytical model, 
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the effective conductivity was assumed to equal the conducti- 
Weary, ethereby treating it as pure conduction and 


a nan “CON 


Equation (1) was non-dimensionalized using the following 


relationships 


- T 
6 = See (a non-dimensional temperature) 
ea ys 
M A 
t= = (a non-dimensional distance) 
oO 
Eo give 
ao 2 


1 2 
E Og a) oe 
with boundary conditions 


es) 


tla = 0 at E =. 0 
: de _ Bs Bes 
and ae =-8(6 - sin wt)at € = 1 
hr 
where B= = is the Biot modulus (which compares the 


pelative magnitudes of the effective heat transfer coeffi- 
erent across the air gap and the internal conduction 
Fesistances to heat transfer). 

An initial” condiiion was not specified as the only 
concern was with the steady state, periodic behavior. Fol- 
lowing Arpaci [Ref. 1], a complex temperature was defined as 

* . 

y(r,t) Sica, t) +. 20 Or, =) 


where (r,t) satisfied 


Gao ee 
ad E ae: C) 
a aes ie Do a 3) 
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watth boundary conditions 
ty = ae ae 2 = 0 at € = 0 
and ay = Se(ieen y= Sete ces OED 2i-—B(6—sin wt) }at € = 1 
fias leads to @(r,t) which satisfied 
28 ) 2 


i aCe ie “0 a0 
Soe. & mOe 


with boundary conditions 


a = 0 gt = 10 
and “ EioCO — sim wer) at & = 1 


ra ; 
arso 0 (r,t) which satisfied 


with boundary conditions 


de 8 : 

OE = at € = 0 

de” I * 2 
and OE = =6CO —elos: Wt) at —& = 1 


a solution of tthe. form 
Waste) I=) diet Ce) 


was assumed, where for large values of time T(t) was 


assumed to equal aoe. therefore, 
W(r,t) = o(r)e7"* (4) 
Equation (4) was then substituted into equation (3) 
dceddy iwr 76 
i oi aa ee, (5) 
c ae a 
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auola thease r 


> 
i} 
+ 
. 
& 7 * jt § 
5 iw 
sup 
aa/ 
“4 
7 | 
Teeth 
ik 1h 
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c 


with boundary conditions 


se = 0 | eee 
Bait a Pep (jOd) Sere 


Equation (5) was expanded to give | 


Bi perem ad oes. 
de? is E ae a 
iwr ? 
Now, let Z = a 3 


and substitute into equation (6) 


af = 0 ae Zz 


and 


i] 
(o) 


| 
a 
: 


The general solution of equation (7) is 


d = Cy I (4) oe Cy Ky (Z) 


as given in Ref. 13 with 


a iG 
(52) 
Io (2) = 14(g2)*4+ 4-5 + ---- = 
s (2!) 
and 2 
Rt4)) = ~{y+log($Z) }1, (2) + >» 
n=l 


Now using the boundary condition 


dd 
Se = 0 at Z = 0 


vi 


(6) 


(7) 


(8) 


ay 


=~ ~Ay A ¢ 
~- ba ! > o°b 


SSbpaq¢e daw. Coy em 


j 


Jo Ast 


arotstibeos waegoe 


om —— 7 = 
+ 
yt} 
1} 
xz t Fits ‘O33 
» @ 
| $4 
? 
> 
| . *Y ; + 


and differentiating equation (8) yields 
‘ d (15 (Z)) d(K,(4)) 
a er, eee 


where 

d(1 5 (2)) 

7 an 0 at Z = 0 
and 

A(K(Z)) 

a) # 0 | atvzZ =. 0 
therefore C, = 0 
and od = Cy I (4) 


Now using the second boundary condition that 
= 


iwr 
(IEE agente rd Santee 
qZ 5 Cd=d) - at Z a 
iwr 
oO 
a 
and differentiating equation (9) gives 
@ GiaG2)) 
do => 1G g 
Z 1 QZ 
2(12(2)) 
Noting that Gann a ys I, (4) and substituting into 


equation (10) 
ae 2 
cI A; Capea © Sa 
et O ) oes 
Os, 
a 


Rearranging and solving for C 


i 
Cy = 1 
iwr | iwr * iwr ? 
o fo) nud fe) 
ag? a a 
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(9) 


(10) 


C¢s) 
gee 
iy en 


@ud then substituting into equation (9) 


~ (2) 
d = — Cid) 
(yee Le tre = 
be a ) i: B O, 2) \ ro} 

Now as 

J, (imx) = i 1 Gi) hRet 3), ps5] 

wr ? wr 2 

- (y ° ie) 0 (J a5 3/2 , 

oO a Oo a. 
and 


wr ? ( wr * 
Or) fe aie y S372 
Ey. ( a 5 :) ro sk oO - :) 


wr * 
Let a = \ — and substitute into equation (11) 


Q 
Jo4 ie af) 
od = ea See (12) 
rene = 372 mc a) 
Seoul 
mee = oe! G = cos 3 + 4 sin an eee (he 
v2. 
substituting this into equation (12) 
eta af) 
ee C13:) 
arene Lames ate Se (144)5,¢0'" =) 
(0) Y28 


ee Wl 
As TaCagi*/ 4) J faze” 4) 


BER (aé)+ iBEi (aé) 


and 


Rey) 
Ce a 


aye J, (abe* 4) 


i 


BER, (a&)+ iBEi, (aé) 


Substituting these results into equation (13) yields 


BER Gacy + LIBEL Caé) 
Ps ce) ° $$ ——— (14) 


. 2 a 
BER (a) + iBEi (a)+ Vag (1-4) (BER, (a) + iBEi, (a)) 
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\e 
ic a* 
aoe ae 
P a 
G 
rf. a x 
> 2( 
gf 
-—_ ap>\t 
7 _ t hae 
(flues? sasdy eukiog 
toa r (3a). Age 
= —_ en o «= 
tags “te) 988 


After rearrangement 


BER (aé) + 2B ee 


a = 3 
[BER_(a)+ —=-BER, (a)+ == Bet, (a)|+ {[BEi_(a)+ 4 BEL, (a)- BER, (a)] 
a 28 v2 a V28 28 
Letting | 
X= BER (a) =—=—BER. (a)+ —*—BE1. (a) 
R re) /28 al /28 af} 
and 
ae a 


BEi (a)+ Bei (al) 2a BER Ca) 
fo) B 1 /2B He 


and substituting into equation (14) gives 
BER. (a&) + iBEi_(aé) 
d = (AC SS oe 
Xp, oT ix, 


Rationalizing the denominator yields 


BER, (a) + iBEi (aé) 


Qe CUCXL, OC- O1X~.) (15) 
xX 2 +X 2 R at 
R 
Now 
_ (BER (ab) X,+BEi (ab) X, ti (BEi, (aé)X,-BER (a&)X,) 
x Z 4. ac 
R al: 


which after rearrangement gives 


x eee H BEA,” (ab) is” 


C1G)) 
2 2. 
X, + X 
where 
oa ey BEi (a&)X, - BER (a&)X, 
= tan 


BER | (a€)X, + BEG) (a&)x, 


80 


‘ane 
JAomag: 
i 


ae 
23m 


iG 
hb» 


» 


Substituting into equation (4) gives 


/ Zz See . 
‘been BER) (a&)+ BEA (aé) jilut Aes 
2 2 


Xx, a Xs 


which also equals 


Z 2 
Xp + x, 


: BER “(aby#BEL.* (ak) : : 
W(r,t) = BY — a> ——__ [cos(wt+8 )+i sin(wt+é )] 


As this problem was modeled as a sine wave, the imaginary 
part’ of (r,t) was used. 
2 2 : 

BER, (a&)+BEi |” (aé) * 
Piycr,t))=0'(r,t)= : 2 "9 ——— sin (wt+6d ) (17) 
which is the analytical solution of infinitely long concen- 
Beve Cylinders experiencing a periodic sinusoidal temperature 
Variation on its outermost surface when heat conduction is 
assumed to be radial only. 


In summary 


SR re 
T-T BER (aé)+BEi -(aé) 
A re) oO P % 

O(r,t)=, ae ieee, ei Ae sin (wt+ 6 ) 

Oe can D3 See co 

R al! 
where 
f-— the temperature of a point r in the rocket motor at 
Sy tele? a 
TAF average bulk temperature of the storage container 
Tyo maximum bulk temperature of the storage container 
W = frequency of the sinusoidal variation(27/24 hours) 
t = time 
r ; é ° 
E = = = dimensionless distance from the center of the 
fo) rocket motor 
t= distance to the surface of the rocket motor 
ry = distance from the center of the rocket motor 
Or 

a= rm = conduction parameter 
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= SS = thermal diffusivity 
Cc 


p = density 


k = thermal conductivity 
© = specific heat 
BER = real Bessel Function 
BEi = imaginary Bessel Function 
a - a F 
Xp = BER (a) + BER, (a) + BEA, (a) 
V28 /2B 
P a ; a 
xX, = BEG (a) +f BEi, (a) - — BER, (a) 
hr | ¥28 V26 
B = cs Biot modulus 
; = BEi (a&)X, - BER (ab)X, 
o => tan 


BER (ad) xX, + BEI (a€)X, 


The following computer programs were used to investigate 
@ wide variety of parameters. The outputs are samples of 


some of these parameter studies. 
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APPENDIX C 


TRUMP Solution 


TRUMP is a computer program for solving transient and 
steady-state temperature distributions in multidimensional 
systems. This program was developed in 1965 at the Lawrence 
Bagration Laboratory by “As L. Edwards [Ref. 2] for their 
CDC/3600 computer. The program was adapted to the Naval 
Postgraduate School IBM/360 Model 67 computer system in 
971 “by Cl. Exbayrum [Ref. 3} from a version used by the 
Be PF. Goodrich Corporation. 

TRUMP is a multipurpose program able to solve a wide 
Varrety of problems involving flow in various kinds of 
potential fields such as heat flow in a temperature field. 
TRUMP allows the solution of general nonlinear parabolic 
partial differential equations both in steady-state and 
transient problems. Complex geometric configurations with 
multidimensional flow may be solved using various coordinate 
Systems. Initial conditions may vary with spatial position. 
Material properties, boundary conditions, and other problem 
Parameters may vary with spatial position, time, or the 
primary dependent variable. | 

ipwe wdatalane fed to TRUMP in “Block" ferm through its 
2 input data blocks, <A complete description of each of 
these blocks is piven in Ref. 2.. A model of the problem 
must be constructed and data: from this model read into 


TRUMP through the data blocks. 
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Two models were used to simulate the rocket motor stor- 
age container system and several variations of each model 
were investigated. 

The first model assumed one dimensional heat transfer 
(radial) with the assumptions that the system was infinitely 
long and that the container surface temperature was spatial- 
ly uniform. The system was modeled as two infinitely long 
Eem@eentric cylinders separated by a 2.94 inch air gap. The 
inner cylinder was constructed of 4130 steel and was filled 
With dry wind blown sand. The thermal properties of the 
materials used in the experimental system are given in 
Table II with units most easily compared to the actual data 
obtained from the system at China Lake. 


TABLE 2 
Thermal Properties of Materials 


Material Density Speciiic- Heat Thermal Conductdavity 
Sand 0.054861lbm/in  O.195BTU/1bm°F 0.00026BTU/min-in-°F 
Steel 0.2807 1lbm/in? 0.109BTU/1bm°F 0.364 BTU/min-in-°F 
Air 0.0000436 " 0.240BTU/1bm°F 0.0000225 Y 


The model was subdivided into volume elements or nodes 
with the representative nodal points given in Figure 27. 
Although the representative nodal point may be located any- 
were in the mode or on the surface of the node, in transient 
problems it is usually located so that the lines connecting 
the nodal points are perpendicularly bisected by the con- 
nected area. This gives maximum accuracy. Two boundary 
conditions were given to the surface node. The first was a 


Sinusoidal disturbance which closely modeled the actual 
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Foeeure 27. “Location of Nodes for One Dimensional TRUMP 
Model. 


Oa 


average experimental data obtained at China Lake as given in 
Appendix D. The wine wave exhibited a maximum temperature 
@eei36°-F and’ an average temperature of 1O4°F. Its period 
was 24 hours (1440 minutes). The second boundary condition 
was the actual average surface temperature of the storage 
container given at two hour intervals. Both these boundary 
conditions are approximations of the actual surface tempera- 
ture. Two hour intervals were the minimum allowable for the 
tabular data as this version of TRUMP has a maximum table 
See. Of 12. 

Several assumptions were initially made. The thermo- 
couple data obtained from the experiment at China Lake gave 
awenaverage temperature at a’ point on the storage container 
and not the actual outside surface temperature. As this 
container wall was only 1/16 of an inch thick and made of a 
good thermal conductor, it was decided to model this data as 
a zero volume boundary node with a known temperature impressed 
on it. It was also assumed” that heat’ transfer across the 
air gap occurred only by radiation and conduction, neglecting 
the secs of free convection, 

It was estimated that the surface emissivities for the 
recket motor and) the’ storage! container were 0.9 [Ref. 6] 
based on their haze gray surfaces. The radiation exchange 


factor for this geometry is given by 


eo, as 
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A sample input deck for the tabular approximation of the 
boundary condition is given at the end of this appendix. 
Several cycles of output data for the one dimensional model 
are also given. 

The second model assumed two dimensional heat transfer 
(radial and circumferential) with the assumption that the 
system was infinitely long. The same physical model was 
assumed for the system except 48 nodes were used instead of 
12. The representative nodal points are given in Figure 28. 
Tae four surface nodes (12, 24, 36, and 48) each had two 
different temperature approximations applied, a sinusoidal 
representation and a tabular input taken at two hour inter- 
vals. The four surface nodes were also modeled as zero 
volume boundary nodes. Each internal thermal connection 
between nodes is described in the input data by specifying 
the two node identification numbers, two connector lengths, 
and two interface dimensional factors. An example of the 
thermal connections ae node 4 is shown in Figure 28 and the 
input data in BLOCK 5 of the two dimensional TRUMP program. 

imencaLculatiton of the radlosities in the two dimension- 
al case was accomplished by using a radiation-network and 
the mocita of crossed-strings. sl 

The radiation shape factors. for the two dimensional sys- 
tem were determined by the method of crossed-strings [Ref. 
ji wine eraphical construction for this method is given 


in -Fieure 29). 
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Pagurce 26: Location of Nodes for Two Dimensional TRUMP 
Model. e 


Figure 29: Graphical Construction for Crossed-Strings 
Method . 
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Be Eo is defined as the fraction of energy leaving sur- 
face m which directly reaches surface n. From the physical 
dimensions of the model A)= A,= A,= 9.45 Bese A,= A,= A, 
= 14.05 $9252, Ay = 12.65 de En, assuming unit depth. Let 
S equal the Length ‘of Aj. 

The crossed-string method lets each surface represent 
the effective surface obtained by stretching a string tight- 
ly over the radiating face between the bounding edges, to 
produce a surface that cannot see any of itself. For ex- 


Bumvey surface 2" in) Figure 29 stretched over surface 2. 


it, whieh by reciprocity leads to 


M1 


By definition F 


22 
oy emer pip = O-2 
2 
and therefore since 
Foo + Foot = 1] then Foo = "0.10 


To calculate the direct radiant heat exchanged between 
surfaces 1 and 2, a minimum-length line was stretched con- 


mecting edge B of A, to edge E of A, and a second ninimum 


2 
to edge F of A 


1 


Henoth line irom edge... -of “A These lines 


dL a 


are labeled Ly in Figure 29 and are equal to the width of 


the air gap, Ly = 2.9375 in. Minimum length Lines were also 


Stretched from point B on Ay EO F on A, and L on Ay fo: E on 


A The length of these lines is Di and is made up of two 


he 
parts; Di'> the tangential distance from F to surface Ay 
and Dis the arc length from the point the tangent hits Ay 


= cay 7h Se " 
to B. From geometry Dy vr tT) 6.62 


Dee sr Week = Ga G2s: 
a 
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therefore 


= ‘ Lt" W 
Di Dy ae Dy Ld 6 O% 
Dae 26 
= 1 Lies 
Now Fio = 7S = 0.86 
1 
Prom reciprocity, Ay Pio = Ay Po4 
12 
Riso SF eno, 528 
Dit x 12 ; 
eli 
Now Fi3 Hse cavlecudated nom 
Fio + 2F,3 = ] 
therefore Fi3 = 0.07 
From symmetry no ae Fi3 


apd then by xeciprocity 


18 
ee. = 
Poy ots Fi 0.047 
il 
Now F was calculated by stretching minimum length lines 


23 


Eton to EL, from E to"*G, from F to/C and.from E to E; where 


the eleneth of the, Lime from oF ‘te. EL= from E to G = S,', 


z. 
from F te G = 2D), and from, Ete, E = 0 
DS a yee 2) 
2 ui 
2 
Fos was now calculated from 
ope tog T Roy t E55. 


therefore 0.002 


Pos 


As Fos was much smaller than the other Fae it was not 
included in the radiation-network diagram in Figure 30 


pRef< 5} 
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Figure 30: Radiation Network. 
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HoOmcabewhate: F 


on the blackbody 


n’ 


set to unity and all other blackbody 


zero. An energy balance was written 
set of six simultaneous equations as 
Node 1 
A, Q-J,) = ALFy (5 ,-J,)t 2A,F 13 
Node 2 
A, (0-J,) = ALFin (Jar -J pt 2A, Foy 
Node 3 
1g) US ES el et Camere icy 6 
ALF, 36537 -J 6+ A, Fi30J 
Node 4 
ee ag) AF 5 y 
Node 5 
A, (0-J,) = 2A,Fo4(Je-Ja)t2A,Fo, 
Node 6 
A, (0-J_) 2A,F13656> «J aus ALFio 
where J. = radiosity of node n 
= 
baa 
eu ie oA 
al 
& 
2 
eee ee Fo 
2 
Lae = radiapion shape factors pr 
Now the values of A., A, and F 
a} 2 nm 


energy balance equations which were t 


fOLM asmshnown in Table Lil. 
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potential of area 1 is 
potentials are set as 


at each node giving a 


follows: 

(J,-J4) 

(J,-J,)+ 2A,F,,(J,-J, 
375: 9+ A,F 136337J pt 
37 5" 

J,-J,)+ AoFo,(J,-Js5) 

oe te gig? 

(J¢-J«) 


eviously calculated. 


were substituted into the 


hen put into matrix 
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betting €,= € 


1 Dee 0.9, a standard computer solution for matrix 


problems gave the radiosities as lictred t2nlabile IV. 


TABLE LV 


Radiosities at Nodes 


{i 


-9046 
~0529 
-00495 
-000797 
~-000125 
-000080 
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Now jo find the radiation exchange factors from node 1 to 
modes n, the radiation network shown in Figure 30 was re- 
duced to the equivalent network shown in Figure 31. Where 
the nodal equations are 


S 


A, (1-0) = E,A, (J5-0) 


1-2 
where eae 
J = eae Oe =e eel: 
1-2 Ay 2 
1341 (1-9) = E,A, (34-9) 
lhe JN 
RG 
Ee ee mn 


These values of the radiation exchange factor are used 
in the two-dimensional program. A sample input deck for 
the sinusoidal boundary condition is included at the end of 
ties appendix. Several cycles of output data for the tabu- 


lar boundary condition are also given. 
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Figure 31: Equivalent Radiation Network 
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APPENDIX D 


Experimental Data 


The data presented in this appendix were obtained from 
the thermocouples on the rocket motor storage container sys= 
tem located at China Lake, California. The thermocouple 
Output was read out on a Honeywell Electronik 25, 24 channel 
Becorder which had been calibrated at 50, 100 and 150°F. 

The data was taken on two consecutive, typical summer days 
fenmemst L and 2, 1972) at China Lake. Each thermocouple 

was read once every 24 minutes. The first set of data pre- 
Senta Ehe storage container temperature at four locations 
plus three different ways of averaging this data. It also 
presents the ambient temperature and the approximate time of 
day. The second set of data presents the surface temperature 
of the rocket motor and three ways to average this data. It 
also presents the temperature at the center of the rocket 
motor and the approximate time of day. Figure 32 shows the 
location of the thermocouples used to collect this tempera- 
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Figure’ 32: Thermocouple Locations for Experimental Data. 
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APPENDIX E 


Uncertainty Analysis 


ag uncertainty analysis was carried out om both the ana- 
lytical solution and on a one dimensional TRUMP model of the 
rocket motor storage container system. In both models, the 
Zmoumetric heat capacity of the sand (9), the conductivity 
of the sand (k), and the emissivity of the surfaces were each 
varied by ten percent to determine the sensitivity of the 
System temperature response to each variation. Although 
other factors may also be varied, it was theorized that these 
tixree had the greatest effect on the heat transfer of the 
system. These factors were also known with the least accu- 
racy; the maximum uncertainty of each was estimated to be 
plus or minus ten percent (odds 20 to 1). 

In the analytical solution, varying the volumetric heat 
capacity changed parameter a, varying the emissivity changed 
parameter 8, and varying the conductivity changed both para- 
meters a and §. The effects on each parameter from each 


Variation are given in Table V. 


TABLE V 
Change in Parameters due to Changes in Thermal Properties 

Change in Property Change in Parameters 
Volumetric Heat Capacity + 102 Eg ae 

Volumetric Heat Capacity - 10% Bye lee 

Emissivity + 102 B+ .49 

Emissivity - 10% ee ccmrnepe 7 

Conductivity + 102 ave ele |Prua 5 22 
Conductivity - 102 Se oe oe 6 oo 
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Each factor was varied holding the other factors constant. 
The changes in temperature and time delay were computed from 
the difference between these new values and those previously 
obtained from the analytical.solution. To obtain uncertainty 
bounds on the analytical curve, the second power equation 


[Ref. 16] was used, namely 


OW. = resulting uncertainty in the calculated temperature 


due to uncertainties in temperature caused by 


Wo = estimated uncertainty in volumetric heat capacity 
= estimated uncertainty in conductivity 
W = estimated uncertainty in emissivity 


naeecoentical calculation was carried out to calculate the 
uncertainty in time delay. The results of these calculations 
are wsmwown in Figures 12 and 13 for the surface and center of 
the rocket motor respectively. The uncertainty in tempera- 
ture varied with time with a maximum variation of +2.75°F at 
Ehe center of the motor and. a maximum variation of +1.85°F at 
the Surface of the rocket motor. The time delay varied by 
oe minutes at.the, center ef the motor and +11 minutes at the 
surface. The actual experimental data was also plotted on 
these Figures for comparison. 

The experimental data also had an uncertainty bound. 
ihree primary factors made up this uncertainty bound; the 
accuracy of the thermocouple wire (+1.5°F), the readability 


oe the recorder (t1°F), and the variation in temperature 
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caused by inaccuracy in the placement of the thermocouples 
ci’r, estimated). The overall uncertainty in the experi- 
mental data was also calculated from the second power equa- 


felon 2s 


~ 


= f 2 2 2 
On SiR REDS PLACE 


DME 
Miese uncertainty bounds are also shown in Figures 12 and 13. 

eB procedure, similar to that used to find the uncertain-— 
Etes,ot the analytical solution, was used to analyze the 
Besulting uncertainty in the TRUMP numerical calculation, 
The results of these calculations are shown in Figures 14 
and 15. The uncertainty in temperature varied with time 
with a maximum variation of +2.95°F at the center of the 
rocket motor and a maximum variation of +1.95°F at the sur- 
Eacenot the motor. ihe time delay varied from +20 minutes 
Jeeenecenter of the motor to +9 minutes at the surface of 
the motor. 

On the basis of the propagation of uncertainty analysis, 
it was determined that the solutions were most sensitive, in 


Onder jof importance, to changes in the volumetric heat 


Capacity, emissivity, and the conductivity. 
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ar to remain stable under desert type conditions and produce 
liantly colored displays of the temperature field. 


a —EEE 
my 314/73 (Wace) 
1 -807-6811 167 Security Classification 


4-31408 


zz jtotsue sleek meet 
i iijolek @ os ee LT: , 


byt Yinvetgaly 


| | 
) 
f 


careht ys tom ie ‘ty | eh 


_ 7G o. 4a| al har be Pie Gee ae cnet 
. = wal z 
iéy@ 6A & ~” 7 onee ” toslenthine ae 


batiwilew @ologidbyserh -eiselay 
a a “inten? as) 
yt ,9 Is ’ | 
rotttled rots if 
<olau Jed vw he « jclisteerkd> elena 
bine Laetoyleda gaeleu be sea) )0Sver eee 
‘bea j ray ( sieta A100 oCidautiagg 
; ’ ( vat ieil’eioe wages age 
Peoy tA a ispig @ebeareqaes 
»\ iw iy aj¢9 SAR 40e@ li ezveate 
; eau Tadnegitages See 
ie sv. 9( 05 Gdo0aeios Jandaieeee 
io bshevé8- daw Wmbgesgloseune 
‘yco bheptt sivetesiasa 
7, | ‘els08 Bee GER (itow 
. Liban / 1146608 vobwi atdags 
ae wet of? So eypeleuli Baga 


Security Clarsificetion 
TALE SRE SERS BE! NMS SR FEI I ST LEE ER ALIS STS Pe EE TEE 5 SINTERS PCNA LE IOC IE AS TL PIE LID MCLE COLE LIS SEITE PRN 


LINK A LINK B LINK C 
KEY WORDS ‘| 


ROLE ROLE WT ROLE 


wr 


meat transfer 
conduction 
radiation 
convection 
peomeentric cylinders 
TRUMP 

Pruawad crystals 

dump storage 


environmental effects 


eis IEE verdes PSA ALES aad Pane a SS 


1 Nov 21473 (BACK } — 
01-807-6821 Security Classification A-31409 


168 


a 


8 Sin 


” 
sna 


Sa 
comer 


a Ce enemies naaa’ edie ne 


‘ 
a Rn se 


Sirs Se 


; Pap fh? ys = a JA it 
See 
ar an et)? 


sio0tls [ndnemmrgsty 
= tte: aan ie 


Se 


141853 


Wirzburger 
An environmental heat 
transfer study of a 
rocket motro storage 
container system, 


thesW6525 
Se ee je : 


An environmental hea 


/ Ill i | I iil te 


“DUDLEY KNOX LIBRARY 


eo 


